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I. INTRODUCTION 

Experimental and theoretical investigations of the 
spin-orbit coupling in semiconductor nanostructures have 
attached considerable attention. 

The Rashba spin-orbit interaction in law-dimensional 
semiconductor structures has been the object of many 
investigations in recent years. The Rashba interaction is 
of the form [Eg 



Vr 



TR\0xPy ~ VyPx) 



(1) 



with standard Pauli spin-matrices o~ x and a y . The 
Rashba interaction can be strong in semiconductor het- 
erostructures and its strength can be controlled by an 
external electric field. 

Recently, advanced growth techniques have made it 
possible to fabricate high quality semiconductor rings. 
The circular quantum rings (nanorings) in semiconduc- 
tors can be described as effectively two-dimensional sys- 
tems in confining potential V c (p) (p = y/ 1 x 2 + y 2 ). The 
parabolic and infinite hard wall confining potentials are 
the most commonly used approximations for the quan- 
tum rings 0. As it was noticed in [ij], these models do 
not allow us to study the tunneling effects and do notper- 
mit the existence of unbounded states. In papers (HQ; 
a more realistic model which corresponds to a quantum 
ring with a potential well of finite depth V = constant 
was proposed: 



V, < p< p t , 

V c(p) = { 0, pi < p< p , 
V, p < p < CO, 



(2) 



where pi and p are the inner and outer radii of the ring, 
respectively. However, this potential was not used in the 
presence of the Rashba spin-orbit interaction. For exam- 
ple, in paper Q, the Rashba interaction was considered 
for the infinite hard wall. In this model V = oo. 

In the present paper, we examine a model which cor- 
responds to the Rashba interaction (1) and the confining 
potential (2). We shall present the exact coordinate wave 
functions and energy levels for this model. 
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II. ANALYTICAL SOLUTIONS OF THE 
SCHRODINGER EQUATION 

The single-electron wave functions satisfy the 
Schrodinger equation 



(3) 



^ + V c {p) + Vr ) * = EV, 



where p is the effective electron mass. 

The Schrodinger equation is considered in the cylin- 
drical coordinates x — p cos tp, y — p sin ip. Further it is 
convenient to employ dimensionless quantities 



2/i 2 P 2 P 2 T/ a 



2/t 



r, 2 

P_ 

P» 



r, 2 ' 



PoPr, 



Pi_ 

Po' 



(4) 



Note that 1 — r» = (p — Pi)/ Po is a relative width of a 
ring. 

As it was shown in 0,9, Eq. (3) permits the separa- 
tion of variables 



* m (r,v) = u{r)e m * ( + w(r)e 



m = 0,±1,±2, 



due to conservation of the total angular momentum L z + 

In the examined model, we look for the radial wave 
functions u(r) and w(r) regular at the origin r = and 
decreasing at infinity r —¥ oo. 

We consider three regions < r < ri (region I), Ti < 
r < 1 (region 2) and I < r < oo (region 3) separately. 

In the regions I and 3 (v > 0) we have the following 
radial equations 



r —r + r— (kjk r + m )u 

ar 



, d w 



dr 2 dr 



dr 
dw 



dw 
dr 



m + 1 



{k+k a r 2 + (m + l) 2 )w 



(6) 



2 



where 



where 7 is defined as follows 



P 2 , .0 

K (e,v,P) = \jv - e- — ±i~. 



(7) 



Note that in the case (3 — (fc+ = k~) Eqs. (6) are 
the modified Bessel equations. In the regions 1 and 2 we 
must select the different particular solutions of Eqs. (6) 
in order to reproduce the correct behavior of the wave 
functions at r — ¥ and at r —¥ 00. 

In the region 1 using the known properties Q 



■%- + -) In(kr) = kln-iQer), 
ar r 



I n {kr) = kl n+ \(kr) 



(8) 



of the modified Bessel functions I n (kr) it is easily to get 
the exact solutions 

ui{r) = cifi(m,r) + digi(m,r), 

wi(r) — —cigx{m + 1, r) + di/i(m + l,r) (9) 

of system (6) where 

fi(m, r) = ^ (I m (k-r) + I m {k+rj) , 
gi (m,r) = l - (l m (k~r) - I m {k+r)) (10) 

are the real linear combinations of the modified Bessel 
functions with complex arguments. Here c\ and d\ are 
arbitrary coefficients. The radial wave functions Ui(r) 
and w\ (r) have the desirable behavior at the origin. 
In the region 3 using the known properties [9j 

d n\ 

— + - K n (kr) = -kKn-^kr), 
ar r J 

^--\K n {kr) = -kK n+l {kr) (11) 
ar r J 

of the modified Bessel functions K n {kr) it is simply to 
derive the exact solutions 

u 3 {r) = c 3 f 3 (m,r) + d 3 g 3 (m,r), 

w 3 (r) = c 3 g 3 (m + l,r) - d 3 f 3 (m + l,r) (12) 

of system (6) where 

f 3 {m, r) = ~ (K m (k-r) + K m (k+r)) , 

g 3 (m, r)= l - (K m {k~r) - K m (k+r)) . (13) 

Here c 3 and d 3 are arbitrary coefficients. At large values 
of r the functions /a (to, r) and .93 (to, r) behave as 



h{m,r) 
g 3 (m,r) 



7T e -W»- e -P 2 /i f[3r + 



■ cos 



7 



2 (v-e) 1 / 4 ^ V 2 
2 (v-e) 1 /^ sm 



,(14) 



y/v-e-P/4 . (3 
cos( 7 ) = == , sm( 7 ) = (15) 



27^6 



Thus, the radial wave functions u 3 (r) and w 3 (r) have 
the appropriate behavior at infinity. It should be noted 
that the radial wave functions have the infinite number 
of zeros for the finite value of energy. 

In the region 2 (v = 0) the radial equations may be 
written in the suitable form 

2 d U dll 9 2\ 

r —5- + r— + (k J k, r — to )u 
ar 1 dr 



dr r 



w 



,dw dw 



r~-r^- + r^f- + {kfk: r 2 - (to + l) 2 )w 
dr 2 dr 

du TO 



< k t-KV[ Tr --u 



where 



(16) 



(17) 



Now in the case (3 = (kf = k i ) Eqs. (16) are the 
Bessel equations. Therefore we use the known properties 



9] 



~r + - I Jn(kr) = kJ n -x(kr), 
ar r 

£ - 2 ) J n {kr) = -kJ n+1 {kr), (18) 
ar r 



^ + -)Y n (kr) = fcy n _i(*r), 



y n (fcr) - -fcy n+ i(fcr) (19) 



of the Bessel functions in order to obtain the exact solu- 
tions 

"2(f) = C21/21 (m,r) + ^21521 (to, r) 

+ C22/22(m, r) + d 2 2ff22(TO,r), 
W2(f) = C2l52l(TO + l,r) + d 2 l/2l(TO+l,r) 

+ c 2 2522(to + l,r) + d 22 f22(m+ l,r) (20) 
of system (16) where 

hi{m.r) = - (j m (k^r) + J m (fc+r)) , 

521 (to, r) = i(j m (fcr r )_ j m (fc+r)), (21) 



/ 22 (m,r) - -(Y m (k-r)+Y m (k+r)), 
922(m,r) = \{Y m {k-r)-Y m {k+r)). (22) 
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Here c 2 \,d 2 \, C22 and c?22 are arbitrary coefficients. 
The continuity conditions 

ui(n) = u 2 (n), u'i(rj) = u' 2 (n), 

wi(n) = w 2 {ri), w[(ri) = w 2 (n), 

u 2 (l) = u 3 (l), t4(l) = «3(1). 

w 2 (l) = «*(!), «/ 2 (l)=t«3(l) 



(23) 



for the radial wave functions and their derivatives at the 
boundary points r = ri and r = 1 lead to the algebraic 
equations 



for eight coefficients, where 



X = {ci,di,c 2 i,d 2 i,c 22 ,d 22 ,c 3 ,d 3 ), 



and M(m, e, u, /3) is 8 x 8 matrix 



(25) 



M(m,e,v,P)=^l (26) 



M(m,e,v,P)X = 



(24) 




/i(m,rj) g 1 (m,r i ) 

f[(m,n) g'i(m,n) 

-gi(m + l,n) fi(m + l,ri) 

-g'i(m + l,ri) f[(m + l,ri) 



-f2i(m,ri) -g 2 i{m 1 r i )\ 
-f 2 i(m,n) -g' 2 i(m,ri) 

-g 2 i(m + l,ri) - / 2 i(to + 1, I' 

-g 21 (m + l,r 4 ) -f 21 (m + 1, r;) / 



M 2 = 



- /22 (m, r*i) -522(m,ri) 

-/ 22 ( m : r ») -g' 22 {-m,ri) 

-022(m + l,ri) -/ 2 2(TO + l,rj) 

-.922(™ + 1 ^») -/^(TO + l.Ti) 



/0 / 2 i(to,1) fl2i (m,l) \ 

/^(m,l) ff 2 i(m,l) 

S2i(m + l,l) / 2 i(m + l,l) 

\ , 921 (to + 1,1) /^(m + 1,1)/ 




/22(to, 1) 522(m,l) -/ 3 (m, 1) -33 (to, 1) \ 

/ 22 (to,1) 522 (to,1) -/£(m,l) -3 3 (to,1) 

<? 22 (to + 1,1) / 22 (to + 1,1) - 93 (m + l,l) / 3 (m + l,l) 

522 (m + l,l) / 22 (m + l,l) - 53 (m + l,l) /^(m + 1,1)/ 



Hence, the exact equation for energy spectrum is 

det M(m, e, t>, = 0. (27) 

It should be stressed that in the explored model the num- 
ber of admissible energy levels is finite for the fixed an- 
gular momentum. Note that Eq. (27) is invariant under 
two replacements m — > — (m + 1) and /3 — > —(3. 

If the energy values are found from Eq. (27) 
then it is simply to get the values of coefficients 
Ci, di, C21, d 21 , c 22 , d 22l c 3 and d 3 from Eq. (24) and 
the following normalization condition 

/•OO 

/ (u 2 (r) + w 2 (r)) rdr = 1. 
Jo 



III. NUMERICAL AND GRAPHIC 
ILLUSTRATIONS 

Now we present some numerical and graphic illustra- 
tions in addition to the analytic results. 

Figures 1 and 2 demonstrate the continuous radial 
wave functions for m = 1, e = 17.88591, /3 = 1, r% = 0.2 in 
the case v = 25 and for m — 1, e = 21.4541, /? = 10, r, = 
0.8 in the case w = 100, respectively. Solid lines corre- 
spond to the functions u(r) and dashed lines correspond 
to the functions w(r). 

Tables 1 and 2 show the dependence of energy e on 
the ring width parameter r, and the Rashba parameter 
(3 for two values of angular momentum number m and 
two values of the well depth v. We see that the number 
of energy levels increases if the well depth v increases. 
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TABLE I: Energy levels for v = 25 
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TABLE II: Energy levels for v = 100 
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FIG. 1: Radial wave functions for m = 1, v = 25, / 
Vi = 0.2. Solid line for u(r), dashed line for w(r). 



= 1, 




FIG. 2: Radial wave functions for m = 1, v — 100, /3 = 10, 
n = 0.8. Solid line for u(r), dashed line for w(r). 



IV. CONCLUSION 



In our opinion the examined exactly solvable model 
with the realistic potential well of finite depth is physi- 
cally adequate in order to describe the behavior of elec- 
tron in a semiconductor quantum ring of finite width 
taking into account the Rashba spin-orbit interaction. 
Further we intend to generalize our consideration by in- 
cluding the magnetic field effects. 
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